An oriented process induced by dynamically regulated energy barriers 
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A novel mechanism for the appearance of oriented processes is investigated with a 
flexible dynamical system overcoming barriers. Under non-equilibrium condition with 
external driving, reaction paths deviate from that at equilibrium with an accompanying 
violation of symmetry between the forward and the reverse paths. Although we never 
introduce any external switching of potentials to generate the oriented processes, multi- 
dimensional flexible dynamics promote the oriented processes through this symmetry 
violation. Along the reaction paths, bottleneck points are proposed as a rate-controlling 
factor, which determine the direction-dependent activation energies satisfying Arrhenius- 
like law for the rate constants. In comparison, in stiff systems, the oriented process is 
suggested to appear in different manner from this scenario. 
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I. INTRODUCTION 

In thermal equilibrium, free energy difference is considered as the single driving force to generate 
oriented processes. On the other hand, far from equilibrium, biomolccular processes are often said to 
be "efficiently" facilitated. Then other mechanisms to generate oriented processes may be explored in 
non-equilibrium situation. One possible example may be clearly found in a biomotor showing directional 
stepping motion along a filament [1]. The stepping of the motor is a kind of reaction shifting the binding 
site from one to the neighbors. Since its position does not contribute to the free energy, the directionality 
is not explained by the free energy difference between before and after this reaction. As ATP-hydrolysis 
drives the motor far from equilibrium, the stepping reaction may not be explained by the rules that govern 
equilibrium oriented reaction. Some experimental works [2,3] suggest that the effect of ATP-hydrolysis 
for the stepping reaction is to generate the direction-dependent activation energies. Such dependency 
certainly induces oriented process of the motor, but the question is then to understand how such a 
dependency can appear in non-equilibrium. 

The generation of oriented processes is a problem of general importance in biological systems. Protein 
folding is facilitated far from equilibrium in some chaperonins with ATP-hydrolysis [4]. Although the 
folding processes are originally oriented by free energy difference, the possible facilitation mechanism 
would act additively to accelerate the processes. Quick conformational changes of proteins after excita- 
tions [5] might be facilitated in a similar way. Keeping several possible examples in mind, we want to 
propose mechanisms which lead to oriented processes. In this paper, we never introduce any external 
switching of potentials to generate the oriented processes, but, instead we show that multi-dimensional 
dynamics can promote an oriented process provided the system is far from equilibrium. 

II. MODEL 

Many chemical reactions can be modeled by barrier overcoming process. The shift of a particle in one 
dimensional periodic potential, a standard theoretical frame for particle current or diffusion [6], gives 
a model to consider steady state for barrier overcoming process. Brownian ratchet models [7] possess 
similar structure, and the stepping process of biomotor systems, too. While referring to these preceding 
studies, we introduce a new feature that the periodic structure is not fixed but fluctuates dynamically. 
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Consider a rail on which stable binding sites of the particle are arranged periodically. When the rail 
is modeled as a flexible chain, the particle senses the interaction potential modulated by the temporal 
conformations of the rail, i.e. modulated by the dynamics of the system itself [8,9]. 

Non-equilibrium situation is introduced by two heat baths, one of which with temperature T\ is attached 
to the rail and the other with T 2 to the particle. Although this setup seems a little bit artificial, it has 
the clear advantage that a non-equilibrium steady state of the system is achieved. The resulting system 
may resemble to Feynmann's ratchet heat engine [10], but in the present model the internal degrees of 
freedom of the rail are explicitly considered. 

A schematic picture of the model is shown in Fig.l. Independent variables are the position x^ of the 
particle and those of the lattice sites {x^ for the ID chain, where i is the site index. The particle and each 
lattice site interact through the potential V(x h — Xi), which has asymmetric shape with one minimum 
as shown in Fig.l. The particle senses the summation of these potentials, V(xh, {xi}) = J^i V( x h ~ Xi). 
The time development is described by Langevin equations, 



li = -7±i + V 2 l k BTiUt) - K c ^{ Xl - it) + 



Xi-i+x i+ i\ dV(x h -Xi) 
{Xl 2 } fl£ : 



/ o , Tf, c ,.x dV(x h ,{xi}) 

x h = -jx h + x/2"fkBT 2 S >h (t) , (1) 

dxh 

driven by two heat baths (T 2 > T x ). We adopt sufficiently small values of the temperatures T\ and T 2 
comparing to the potential well depth of V. 7 is a friction coefficient and £ a (t) represents Gaussian 
white noise. We use the units as Boltzmann constant ks = 1. Although Eqs.(l) include inertia terms, 
qualitatively the same results are obtained in overdamped cases without these terms. K c is the spring 
constant for the chain with which each lattice site is connected to neighbors, and also to fixed ground, 
with the natural interval length I. 

Hereafter two typical values of stiffness K c are used for the demonstration: at K c = 0.5, the flexibility 
of the chain greatly affects dynamics of the system (flexible system), while at K c = 5.0 the dynamics of 
the internal degrees of freedom for the chain are considerably suppressed similarly to the cases of certain 
Brownian ratchets (stiff system). 



III. RATE CONSTANTS DEFINED FROM TIME SEQUENCE 

A typical time sequence of the particle position is shown in Fig. 2, in which we can observe each 
jump (stepping process) of the particle. As is seen there, the particle steps to both direction with equal 
probability when T\ = T 2 , whereas oriented motion appears under the existence of temperature difference 
(T 2 > Ti). In the following, 'forward' direction is defined for each value of K c as the direction of the 
oriented motion in T 2 > Ti, and 'backward' direction as the reverse direction. The forward and the 
backward stepping processes correspond to their reverse process each other. 

Because most of the successive stepping events are sufficiently separated in time, each process could be 
regarded as a stochastic one. Then, it is natural to consider the probability of the steppings in a steady 
driven state. The long-term observation of the particle position gives an ensemble of data how long the 
particle has stayed at one lattice site after previous stepping event. Consequently we have a "population" 
of waiting states as a function of waiting time. Time evolution of such population could be analogous to 
chemical kinetics, and a rate constant for the occurrence of a step is introduced as follows. 

Let N a (r) denote the number of events that the particle stays at a lattice site without stepping during 
a time longer than r and Np{r) (Nb(t)) denote the number of events that the particle has stepped 
forward (backward) within r. The rate of the forward (backward) stepping np{r) (kb(t)) is defined as 

^ = k f (t)N,(t), ^ - k b (t)N s (t), (2) 

where dNp/dr (dNs/dr) is the number of the forward (backward) stepping events per unit time. 

In the inset of Fig. 2, Kp(r) and k b (t) are shown for the case of K c = 0.5 for T 2 > Xi. It is found that 
these rates are almost constant for sufficiently large r. For all cases examined with changing Ti, T 2 and 
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K c , converged values for Kf{t) and kb(t) were obtained, which implies the process can be regarded as 
an independent stochastic process for long r. In the following, we let up and Kb, which referred to the 
rate constants for "stepping reaction", denote the constant values of kf(t) and kb(t) for large r. 

From the long time average of the stepping motion, the current and the diffusion coefficient of the 
particle are defined as J = linit-xx, (xh{t) — Xh(0)) /t, D = lim^oo {(xh(t) — Xh(0)) 2 )/2t. For all the 
cases examined numerically, J and D converged to finite values. In thermal equilibrium (T\ = T 2 ), 
the particle simply diffuses along the chain without drifting (J = 0) all over the temperature range, 
while J 7^ in non-equilibrium (T 2 > T\ ) . Interestingly, the direction of the drift in the flexible system 
(K c = 0.5) is opposite to that in the stiff system (K c = 5.0) as shown in Fig. 3. The convergence for D 
means that the particle shows normal diffusion not only in equilibrium but also in non-equilibrium steady 
states. 

If the particle steps stochastically, J and D should be related to the rate constants as follows. 

\J\*1{k f -k b ), D^ P{KF + KB \ (3) 

where I is a lattice interval of the chain. As is seen in Fig. 3, these relations hold well especially in 
low temperature. At high temperature, deviation from these relation occurs since the contribution from 
correlated successive stepping (multiple stepping) in small t is not negligible. Thus, it is confirmed 
that the stepping process is well characterized by the rate constants (kf,k b ) at least within some lower 
temperature range. 

IV. BOTTLENECK POINTS AS A RATE-CONTROLLING FACTOR 

For a simple chemical reaction, it is well known that rate constant k obeys Arrhenius law k ~ 
exp(— Ea/ksT), where E a is called activation energy characterizing the reaction rate. In transition state 
theory (TST), the activation energy is purely determined from the information of the potential energy 
surface along the reaction coordinate without considering temperature effect. The minimum height of 
energy barrier, the difference of potential energy at the saddle point from that at the initial stable point 
along the reaction coordinate, corresponds to the activation energy [11]. 

Far from equilibrium, we need to define which quantity can be used as the temperature to characterize 
the process. Because the particle position xh is the important observable for the stepping process, the 
fluctuations of Xh in the staying state at a site would give the temperature. This is approximately equal 
to kinetic temperature of the particle, Th = {x\), because the stepping event is rather rare and J is 
sufficiently small in the range of our analysis. Here (•) means ensemble and/or long-time average. Under 
non-equilibrium condition (T 2 > 7\), T\ takes values such that T 2 > Th > T\ due to the connection to 
the chain, while Th is equal to T 2 (= T\) under equilibrium condition. 

Let us consider the activation energy for the present system according to TST. From the fixed point 
analysis for Eq.(l) [12], the stable and the saddle points can be determined so that the potential dif- 
ference AUo between them is obtained straightforwardly. Note that AUq is a constant independent of 
temperature, while depends on the flexibility of the system. In the flexible system, there is a cooperative 
effect between the particle and the rail to make the activation energy lower. In the thermal equilibrium 
cases in Fig. 4, the obtained AUq actually works well as the activation energy in the Arrhenius law, char- 
acterizing the stepping rates for the broad range of temperatures. However, once the system is put under 
the non-equilibrium condition, the rate constants deviate from the Arrhenius relation in equilibrium. 

In the following, we try to get a modified version of the Arrhenius relation, applicable to non-equilibrium 
situation. In TST, the saddle point is the most important for the calculation of the activation energy. 
Although the usage of the saddle point itself might not be valid far from equilibrium, the idea that some 
bottleneck point in phase space is the rate-controlling factor should not be necessarily abandoned. How 
can we find such a "bottleneck point" ? The paths for the stepping process is thermally fluctuating, but 
the ensemble of these paths would be recognized as tube-like clouds in a phase space. Then the most 
probable reaction path (MPRP) could be defined by averaging over the tube-like clouds. We suppose the 
bottleneck point would be defined as the saddle-like point along this MPRP. 

In MPRP the particle position Xh, which is the major observable for the stepping process, would be 
the important coordinate in the phase space. According to this idea, we propose that the ensemble 
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which determines the bottleneck point consists of the points in the paths satisfying the conditions for the 
potential derivatives by X] x : 

dV(x h ,{xi}) n d 2 V(x h ,{x t }) 

dx h ~ U ' dx\ <U ' W 

with the supplemental condition that the points have been passed almost synchronously with the actual 
stepping events. The condition (4) is more gentle than that for the saddle points because it neglects the 
derivatives in other directions {xi}. 

The first condition of Eqs. (4) means that force acting on the particle is zero at the point and the 
second condition means that the point is an unstable fixed point if looking only along Xh- By averaging 
the ensemble of the points satisfying the above condition with distinction of the direction of the stepping, 
bottleneck points for the forward and the backward stepping processes are obtained separately. 

Interestingly the bottleneck points for the forward and the backward stepping are often different in non- 
equilibrium (T 2 > Ti), while the two points is always identical in equilibrium (T 2 = Ti). It suggests that 
the symmetry between the forward and the backward processes is violated in the presence of temperature 
difference. 



V. DIRECTION-DEPENDENT ACTIVATION ENERGIES 

The difference of the bottleneck points will lead to the difference of the energy barrier. The barrier 
height AUf (AUb) can be determined for the forward (backward) stepping process from the potential 
difference between the stable and the bottleneck point in similar manner to AUo. It is noted that AUf 
and AUb could depend on the temperatures in contrast to AUo, because the bottleneck point is obtained 
as the ensemble average of the points in the numerically simulated trajectories. 

In Figs. 5, AUf and AUb are displayed as a function of the external driving (temperature difference). 
In equilibrium (Ti = T 2 ), AUf is equal to AUb as directly noticed from the identity of the bottleneck 
points between the two directions. This is consistent with the equality of the rate constants kr = Kl- 
In non-equilibrium, the barrier heights varies from its equilibrium values, and their dependences on the 
external driving are different between the flexible and stiff system as is seen in Figs. 5. 

For the flexible system (K c = 0.5), barrier heights grow with the increase of the external driving, where 
the deviation of AUf from AUq is smaller than that of AUb- As clearly seen in Fig. 6, substituting 
AUf and AUb for the activation energies, the rate constants in non-equilibrium can be reproduced by 
the Arrhenius-like form 

( AU f\ a ( &U B \ 

KF=Aexp {-kBT-J> KS = Aexp l~M:J' (5) 

where A is a temperature-independent prefactor approximately common in the forward and the backward 
stepping. This relation indicates that the bottleneck points are major rate-controlling factors and that 
AUf (AUb) is the substantial activation energy for the forward (backward) stepping process far from 
equilibrium. The knowledge of these activation energies is sufficient to characterize the oriented stepping 
process. 

To be emphasized, the direction-dependent activation energies mean that absolutely different paths are 
selected between the forward and the backward stepping reactions. From the dependency of AUf and 
AUb, the reaction path for the backward stepping process is implied to be more sensitive to the external 
driving than that for the forward. This sensitivity may be related to the property of the cooperative 
dynamics to realize the reaction path but this is open as a future problem. 

For the stiff systems, the fluctuations of the rail is small and the cooperative dynamics between the 
particle and the rail are almost suppressed. As a result, the bottleneck points do not draw apart from 
the saddle point with increasing external driving, and the energy barriers defined by these bottleneck 
points only slightly increase compared to AUo (see Fig. 5(b)). Such slight variation cannot account for 
the deviation of the rate constants in non-equilibrium seen in Fig. 4, so the barrier heights AU f and AUb 
does not explain the relative rate Kf/kb- Further investigation for the stiff systems is not considered 
here, but it can be done in parallel to the Brownian ratchet models which has no flexibility of the chain 
[13,7]. 
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VI. DISCUSSION 



In this Paper, we have discussed mechanisms to generate oriented process. Far from equilibrium in the 
flexible system, the direction-dependent activation energies (DAEs) account for the oriented process well. 
For the appearance of DAEs due to dynamical coupling between components, the symmetry violation 
between forward and its reverse reaction paths is a key feature, and it never occurs in equilibrium. 
Thus oriented processes due to DAEs cannot be explained from the equilibrium properties of the relevant 
systems. Instead, we should develop a direct understanding of non-equilibrium phenomena. The fact that 
Arrhenius-like form could be recovered in non-equilibrium steady state might suggest such new direction 
to find some robust structure to describe the system in non-equilibrium. 

Direction-dependent activation energies are also obtained for the stepping process in Kinesin biomotors 
by analyzing the rate constants for the stepping motion in experimental sequences [2,3]. Comparing to 
the simple diffusion process without ATP-hydrolysis exhibited by Kinesin mutants [14] (proposed as an 
equilibrium system), the directionality in the normal functioning state seems to be generated by the 
suppression of the backward stepping process while leaving the rate of the forward process almost similar 
to that in equilibrium [3] . Interestingly, our flexible system also exhibits a similar tendency under external 
driving. 
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FIG. 1. Profile of the model (refer to the text). Potential form is asymmetric in space, 
V(Ax) — JGjtanh(p(Ax — r))/cosh(dAa;), where interaction is confined mostly to the nearest lattice sites 
(ef _1 = 0.25 with 1 = 1). The degree of asymmetry is fixed at p = 10 and r = 0.3, while coupling strength 
Kh — 0.2. In numerical simulation, the chain consists of 7 or 40 lattice sites with periodic boundary condition, 
while 7 is fixed at 0.1. 

FIG. 2. Typical time sequence of x h for K c = 0.5, Ti = 0.02, T 2 = 0.08 (T h = 0.042) and Ti = T 2 = 0.042. 
(Inset): r vs • kf(t), o « b (t) calculated according to the definition in Eq.(2), which show quick convergence. 

FIG. 3. Current J and diffusion coefficient D for T\ = 0.02 with various values of T2. The guidelines are the 
relationship in Eq.(3). J changes its direction from the flexible to the stiff system, where it is positive in xh for 
K c — 0.5 while negative for K c = 5.0. 

FIG. 4. Arrhenius plots, 1/Th vs kf,kb- In non-equilibrium (neq), we fix Ti at 0.02 while Ti is widely changed, 
i.e., the points in (neq) equal to those in (eq) at 1/Th = 50. The dotted lines are the Arrhenius law exp(— AUo/Th), 
where AC/o ~ 0.15 and 0.22 for K c = 0.5 and 5.0, respectively. Prefactor is a fitted constant. For K c = 0.5, 
«F(eq)= KB(eq)~ KF(neq) while «s(neq)< Ks(eq). For K c = 5.0, both ^^(neq) and fts(neq) are suppressed. 



FIG. 5. The increase of the potential barrier AUf and AUb- The barrier height grows with Th — Ti for the 
flexible system but it does not for the stiff system. Upper (Kc — 0.5): Ti = 0.02 (circle), Ti = 0.01 (rectangle). 
Lower (Kc = 0.5): Ti = 0.02 (circle), Ti = 0.03 (rectangle). At the equilibrium point (T h = Ti) for K c = 0.5, 
AUf = AUb but there are slightly larger than AUo- This deviation might be due to 'saddle point avoidance' 
[15]. 

FIG. 6. Exponential relationship between k f (k s ) and AUf /T h (AUb /T h ) for K c — 0.5 (k B = 1). Arrhe- 
nius-like relationship is recovered with regarding AUb as the activation energy. 
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